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Abstract. We propose a universal decomposition of unitary maps over a ten- 
sorial power of C 2 , introducing the key concept of phase maps, and investigate 
how this decomposition can be used to implement unitary maps directly in the 
measurement-based model for quantum computing. Specifically, we show how 
to extract from such a decomposition a matching entangled graph state (with 
inputs), and a set of measurements angles, when there is one. Next, we check 
whether the obtained graph state verifies a flow condition, which guarantees an 
execution order such that the dependent measurements and corrections of the 
pattern yield deterministic results. Using a graph theoretic characterization of 
flows, we can determine whether a flow can be constructed for a graph state in 
polynomial time. This approach yields an algorithmic procedure which, when it 
succeeds, may produce an efficient pattern for a given unitary. 



1 Introduction 

Measurement-based quantum computing 1 1I2I3I4I has attracted considerable interest 
in quantum computing research during the last few years. Different aspects of this 
model have been investigated, including: characterization of the initial entangled states 
I5I6I7I8I9I10I . description and unification of different models II 1I12I13I14I15I16I17I18I . 
investigation of fault tolerant computing within this framework 1 19 20 21 22 23 2 41251 . 
and proposals for implementations I26I27I28I29I30I31I32I33I34I351 . 

It remains an open question however, whether this model may suggest new tech- 
niques for designing quantum algorithms. This is the question we address in the present 
paper. Specifically, we introduce a methodology for the direct decomposition of a given 
unitary map into a one-way model pattern. 

Previously, one would typically start with an algorithm already implemented in the 
circuit model, replace each gate by a corresponding pattern. To transform the pattern 
to a standard form where all entangling operations are performed first, one could then 
use the by-product method |3| or the more general standardization algorithm 1 16 1. In 
this paper, we propose a direct method that is free from any reference to the circuit 
model. We start with the observation that one-way patterns implicitly define a particu- 
lar decomposition of unitary maps into a preparation map enlarging the input space, a 
diagonal map with unit coefficients, and a restriction map contracting back the space to 
the output space, which we call a phase map decomposition. (Note that this decompo- 
sition does not directly correspond to any physical procedure.) 

Since the one-way model is universal, this alternative decomposition is also univer- 
sal; and there is a straightforward procedure which allows us to determine a phase map 



decomposition for a unitary from a pattern implementing the same unitary. Remarkably, 
one can define a reverse procedure as well, which breaks in two steps. First, given a uni- 
tary map, one enumerates such phase map decompositions by constructing the right set 
of coefficients in the middle diagonal map (Lemma|2]and Algorithm^. Then for each 
such decomposition, one verifies whether there exists a matching one-way pattern. This 
reduces to finding for any phase map decomposition a matching entangled graph state 
(with inputs) and choice of measurements angles (Algorithm [2j. To obtain the execu- 
tion order of the one-way pattern, we look for a "flow" for the underlying entanglement 
graph which describes an adequate set of dependent measurements and corrections so 
that the obtained pattern is deterministic (Theorem[0. Using a graph theoretic charac- 
terization of flows, we can determine whether a flow can be constructed for a graph 
state in polynomial time (Algorithms|5]and|3. 

This constitutes a method for the implementation of unitary maps in the one-way 
model which bypasses completely the circuit picture, and may produce more efficient 
implementations in particular instances. 

2 Preliminaries 

We first briefly review the one-way model and the various notions relevant to this paper. 
2.1 One-way Model 

Computations involve a combination of one-qubit preparations N%, two-qubit entan- 
glement operators Eij :— /\Zij, 1 one-qubit measurements M", and one-qubit Pauli 
corrections X { , Z i (where i, j represent the qubits on which each of these operations 
apply, and a is a parameter in [0, 2ir}). A preparation Ni prepares auxiliary qubit i in 
state |+}». An entangling operator entangles qubits i and j. A destructive mea- 
surement Mf is defined by the pair of orthogonal projections (± | j applied at qubit i, 
where |± Q ) stand for ^j(|0) ± e m |l)). We use the convention that (+ a |, corresponds 
to the outcome 0, and {— a \i corresponds to 1, Since qubits are measured at most once 
in a pattern, we may unambiguously represent the outcome of the measurement done at 
qubit j by Sj . The dependent corrections used to control non-determinism are then X^ 3 
and Z- 3 \ with sj € {0,1}. 

A pattern is defined by the choice of a finite set V of qubits, two possibly overlap- 
ping subsets / and O determining the pattern inputs and outputs, and a finite sequence 
of commands acting on V. We require that no command depends on an outcome not 
yet measured; that no command acts on a qubit already measured; and that a qubit is 
measured (or prepared) if and only if it is not an output (respectively an input). 

Let Hi (Ho) denote the Hilbert space spanned by the inputs (outputs). The execu- 
tion of a pattern consists of performing each command in sequence. If n is the number 
of measurements (i.e. the number of non-output qubits), then this may follow 2" differ- 
ent computational branches. Each branch is associated with a unique binary string s of 
length n, representing the classical outcomes of the measurements along that branch, 

1 The symmetric operator of controlled-.^. 



and a unique branch map A s representing the linear transformation from Hi to Ho 
along that branch. (Note that states resulting from a measurement are not normalized.) 
The branch corresponding to the classical outcomes s = • • • 0, will be called the pos- 
itive branch. 

A pattern is said to be deterministic if all of the branch maps are equal: any given 
unitary can be implemented by a deterministic pattern (in the precise sense of The- 
orem [2 below). A pattern is said to be in standard form if entangling operators Eij 
appear first in its command sequence. Any given pattern can be put in such a form, 
which reveals implicit parallelism in the computation, and suits certain implementa- 
tions well [16|. One can then define an entanglement graph associated to the initial 
preparation and entanglement step, with vertices V and edges Eij. By doing so, the 
computation of the pattern can be understood as the direct manipulation of such mul- 
tipartite entangled states (which are called graph states in the case where there are no 
inputs). 

2.2 Geometries and Flows 

In order to find a pattern for a given phase map decomposition of a unitary operator, 
we will need to search for a matching entanglement graph and write the pattern directly 
from the obtained graph. We may do this using flows 1361 . which we review now. 

A geometry (G, I, O) consists of an undirected graph G together with two subsets 
of nodes / and O, called inputs and outputs. We write V for the set of nodes in G, I c , 
and O c for the complements of / and O in V . We write x ~ y and xy e G when 
x is adjacent to y in G , and Eq := Yixyea ^ x v ^ or tne global entanglement operator 
associated to G. (Note that all E xy commute so the order in this product is irrelevant.) 
To extend a geometry into a pattern, it remains only to decide the angles which are 
used in measuring qubits of O c , and the dependent corrections to be used; thus every 
pattern has an underlying geometry which is obtained by forgetting the measurements 
and corrections. 

We give now a condition on geometries under which one can synthesize a set of 
dependent corrections such that the obtained pattern is deterministic. 

Definition 1 A flow (/, =4) for a geometry (G, I, O) consists of a map f : O c — > I c 
and a partial order =4 over V such that for all x £ O c : 

— (Fi) x ~ f(x); 

— (Fz'z) x 4 f(x); 

— (Fiii) for all y ~ f(%) , we have x =4 y ■ 

The coarsest order for which (/, =4) is a flow will be called the dependency order 
induced by /. 

Theorem 1 (|36|) Suppose the geometry (G, I, O) has flow f, then the pattern: 




where the product follows the dependency order =4 of f, is deterministic, and realizes 
the unitary embedding: 



That is: if a geometry has a flow, we may without loss of generality consider only its 
positive branch. Then measurement-based computing downs to projection-based com- 
puting. This will be handy in the formulation of our procedure. 

3 One-way Patterns as Phase Maps 

With our preliminaries in place, we now turn to the formulation of our decomposition. 
Various operators over H preserve the computational basis, up to phase: for example, 
Pauli maps, Z a (defined as Z a \0) = |0), and Z a \l) = e M |l)), and controlled-Paulis. 
One-qubit measurements also map the standard basis of one space to those of another, 
up to a scalar factor. In particular, one has the following simple equations where \x) is 
an n-qubit computational basis state and j is an index for a qubit in |x): 



We will call a map <P : Hv — > a phase map if it is diagonal in the computational 
basis and has only unit coefficients. The typical example of such a map is AZ. It is 
important to note that the above definition depends on the choice of a basis. 

Following one-way model terminology we also define a preparation map Pi^v ■ 
Hi — > Hv that expands the input space by tensoring auxiliary qubits, 



and a restriction map Ry^o ■ Hv — > Ho that projects the space to the output space: 



It is easy to see that the restriction map is the adjoint of the preparation map. 

As we have seen above, measurement and entangling commands in the one-way 
model define phase maps, and hence from the universality of the model we obtain the 
following decomposition: 




(± a |®J® n (|0)|z» = 2- 1 /2| x ) 
(± a | ® J® n (|l)|x)) = ±2- 1 / 2 e^ a \x) 
AZ y |0)|x) =\0)\x) 
AZ y |l)|x) =Zj\l)\x) 



\x) I ► \x) ® I + •••+) 70 



X) 
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Theorem 2 For all unitary U : Hi 
such that: 



Hq, there exists a phase map <P : Hy 



Hy 



U = Rv^o ^ Pi^v 



Proof. We know that hZ and J a are universal |37|, and also that compositions of pat- 
terns having flows, themselves have a flow |36|. Hence there is a deterministic pattern 
V implementing U, and from Theorem^we then have: 

= 2 |cr|/2 YlitoWr* Eg Pi^v 

= R V ^o n,0o z r a ' E c p^v 

= Rv^o $v^v Pi^v 

where <P — Yii^o Z^ ai IlyeB ^Zij is the phase map corresponding to the entangle- 
ment operations and measurement angles. □ 

One can think of the above theorem as a special kind of diagonalization for unitaries 
where one is allowed to inflate the dimension of the underlying space. This will prove 
to be useful for direct programming in the one-way model. We present first a couple 
of examples, clarifying our first set of definitions and preparing the ground for a direct 
proof of the above theorem which does not invoke universality of the one-way model. 
This then leads to an algorithm for pattern design. The examples already hint at the 
construction behind the direct phase map decomposition algorithm. 



3.1 Example 1 



Consider the unitary map J a : — ► ^{1} which decomposes in the computational 
basis as: 
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This decomposition is obtained from the one-way pattern X^M^En which imple- 
ments J a , and has as positive branch the phase map: 

2 1 /2( +Q | 1 A Z 12 : H {h2} — >H {2} 

Factoring out the restriction operator gives the decomposition: 
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where the left matrix is the restriction Ri : ?in 2} — ► ^{2}- 



Note that this decomposition is not physical, since the phase map doesn't directly 
correspond to any operation. Indeed it is a convention to substitute the physically mean- 
ingful equation (+ a | ® 7(|1)|1)) = 2- 1 /2 e -"*| 1 ) with: 

i? 1 (2- 1 /2 e -^«| 1 )| 1 )) = 2- 1 /2 e -«*| 1 ) 

However, the corresponding one-way pattern is a physical procedure that may be un- 
derstood in terms of the above decomposition. 

The phase map here is Z^ a l\Z\ 2 , and the decomposition above can be rewritten: 

J a = R l {Z^ a AZ 12 )P 2 

As we have seen in the proof above, this reasoning is perfectly general; indeed for any 
pattern, one gets the phase map associated to the positive branch: 

n zr> n ^ 

ieo a ijeE 

and writing U for the corresponding unitary (assuming the pattern does compute a uni- 
tary map): 

u = r c ( n zr- n ^ ) p IC 

3.2 Example 2 

Example 1 uses only one auxiliary qubit, and as such is a special case where the required 
number of auxiliary qubits is equal to the number of inputs. This is of course not always 
the case and the general algorithm for phase map decomposition will take care of this. 
We will present exact bounds on how much one needs to expand the computational 
space to be able to obtain the decomposition; however to realize a decomposition as a 
pattern we will need further restrictions. The following example demonstrate this case. 

The shortest known pattern for the ^"-rotation is X^ 2 Z^ 1 M$M^ a E 12 E 23 with 
positive branch: 

(+\2( + -a\lEi 2 E 23 

which induces the 3-qubit phase map <P\xyz) = (-l) xy+yz e lay \xyz) (indeed a diago- 
nal of units) and corresponds to the following decomposition of Z a : 

R 12 D(l, 1, 1, -1, e ia , e ia , -e ia , e ia ) P 23 

where £'(• ••) is a diagonal matrix. Note that some permutations of the diagonal lead 
to other solutions, and most decompositions won't correspond to a pattern. Also if one 
uses only 1 additional qubit one obtains another decomposition with middle map: 

C(V2,0,0,V2e M ) 

which is not a phase map since coefficients are not units. The natural question is whether 
it is possible to generate all such decompositions without using a pattern, and next find 
a pattern matching one of these decompositions. 



4 Direct Decomposition 



There is no a priori reason why a phase map decomposition could not be obtained with- 
out any reference to the one-way model: indeed, it is possible to give a direct approach. 
Supposing one adds n auxiliary qubits to the input space I, a simple calculation shows 
that for each coefficient u in the the computational basis matrix representation of U, 
there will be 2l l/ l~l°l /2l J l = 2"~I 7 I 'slots' to spread over the diagonal of the phase 
map <P (since U is a unitary \I\ = \0\). Thus, finding a decomposition amounts in this 
case to finding complex numbers a;W such that the following two conditions hold: 

u= Y, xil) C 1 ) 
i<2™-m 

2™/ 2 |a;W| = l (2) 

The first equation says that the restriction map R will sum up all the x^'s to give u, 
while the second one asks for unit diagonal elements (note that the preparation map of 
n auxiliary qubits introduces an overall factor of 2~™/ 2 ). 

Lemma 2 If n > \I\ , Equations\J\and^have joint solutions iff\u\ ^ 2™/ 2 ~l / l . 

Proof. Each complex a;W can be seen as a real plane vector of constant length 2~™/ 2 , 
and all one has to do is to choose their angles in a way that they will globally add up to u. 
If one aligns all x^s with u, the resulting sum is at least as long as u iff \u\ ^ 2™/ 2 ~l / l ; 
thus this inequality is necessary for Equations[2and|2]to have joint solutions. 

If n > \I\ , then 2™~I 7 I 2 , so there are at least two terms . We may pick any 
two of them and rotate them at opposite angles ±9. If 9 reaches ^ before the global 
sum matches u, then the corresponding two x^s contribute nothing, and we pick two 
additional terms to rotate. Clearly, at some stage, for some value of 9 the sum will co- 
incide with that of u. □ 

Due to the unitarity of U, \u\ < 1 : so a safe choice of n is one such that 2 n / 2 ~l 7 l > 
1. Thus, n > 2|/| is always sufficient for a phase map to exist. Another consequence of 
the above lemma is that for any given unitary map U on |7| qubits, unless U is itself a 
phase map and also requires no auxiliary qubits, we have n > \I\: then a lower bound 
on the number of required qubits to implement it as a one-way pattern is 2|7|, if at least 
one coefficient of U is larger than ^ . 

(i) 

For a unitary U , once we have fixed n, an output space O, and a choice of x P q 
satisfying EquationsHandlJlfor the coefficients u pq = (p\U\q), the following algorithm 
will enumerate all possible decompositions: 

Algorithm 1 Input: for sets V, I, and O and n = \I C \: 

— a unitary U on Hi; 

— complex numbers { Xpl } 2 =1 ' ' satisfying Equations^\and\2\for each u pq ; 

— a permutation a over {1, • • • , 2™~' 7 '}. 

Output: diagonal elements {dkk}t=i > such that dkk = V2™ x P q , where: 

— the binary representation of p agrees with that ofk after restriction to O; 



— q=k mod 2l 7 l; 

- i = a([k/2^\). 

The elements {dfcfcljtli are the solution of R<PP — [/where ^ = D{dn, d^, • • • ) ; 
due to the simple structure of matrices R and P we derive the above algorithm. 

Note that obtaining a decomposition is not sufficient for the existence of a pattern. 
We now discuss how to extract an entanglement graph and measurement angles from a 
phase map, if there is one. 

5 Pattern Design 

To determine whether a phase map decomposition R <P P of a unitary U has a corre- 
sponding pattern in the one-way measurement model, one wants a graph Ge over V, 
and angles ay for j £ O c such that 

II Zp II AZ jk 

That means for all x in the ^-computational basis: 

d xx = e - i T.oc<*i*j (_i)E JfceB *jf** (3) 

where d xx is the diagonal coefficient of the phase map corresponding to x basis. Based 
on this observation we propose the following algorithm for the above graph matching 
problem. 

Algorithm 2 Input: A phase map decomposition for U — i.e. the diagonal elements 
{dxx} x =\ from Algorithm^ 

Output: either (i) A graph G on V and a set of angles of measurements {aj}j e o c > or 
(ii) no matching graph exists. 

1. For j € {1, • • • , |O c |}, consider the \V\-bit string Zj that only has a 1 at position 
j, and set a such that e~ laj = d ZjZj . 

2. For all j, k, consider the \V\-bit string Zjk having a 1 only at positions j and k. 
Check whether d ZjkZjk = ± e ~ 2 ( Q i+ a 't) (the angles for the corresponding qubit in 
O is taken to be 0). 

— (i) if YES and the sign is —1, return Ejk as an edge in G. 

— (ii) if NO, no matching graph exists. 

From Equation|3] we can see that a given phase map can be implemented by at most 
one pattern. Once an entanglement graph is obtained, we check whether (G, I, O) has 
a flow, in which case Theorem^obtains a deterministic pattern for U. We address the 
problem of finding a flow for the entanglement graph in the next section. 

Although there are exponentially many elements on the diagonal of the phase map, 
testing for the existence of G will query the middle diagonal map only quadratically in 
|/| + n, to read off the measurement angles and the entanglement graph. This in prac- 
tice could accelerate the detection of bad decompositions before obtaining all diagonal 
elements in the phase map. 



When the procedure fails, because the obtained graph does not have flow or simply 
does not exist, one backtracks by: (1) trying a different decomposition given by Algo- 
rithm^ (2) trying another solution from Lemma|2] (3) revising the choice of outputs, 
and ultimately (4) expanding further the computational space. Without any additional 
constraints, it seems that there are many solutions to be checked. One might be able to 
infer additional constraints to Equations ^andEl from the requirement that there be a 
corresponding entanglement graph, reducing the set of phase maps which we consider. 
How this may be done is, however, an open question. 

6 Finding Flows 

In this section, we outline algorithms which allow us to find flows for a geometry 
(G, /, O) in polynomial time. The proof of correctness of the algorithm employs tech- 
nical results of graph theory, which we will only sketch here: a rigorous development 
may be found in l38l . 

6.1 Graph Theoretic Characterization of Flows 

We will begin by showing how flows can be characterized in the language of graph 
theory, in order to make use of the solved problems in that field. 

First, we present some additional notation related to directed paths. In a directed 
path, we write x — * y to represent an arc between vertices x and y. If C is a collection 
of directed paths, we will say that x — > y is an arc ofC when x — > y is an arc in one of 
the elements of C . 

For a flow (/,=<;), we can easily see that the function / must be injective: if f[x) — 
f(y) , we have y ~ f(x) and x ~ f(y) by (Fi), and x =4 y an d y =4 x by (Fiii), so x = y . 
Therefore, each vertex x S V has at most one vertex w S V such that x = f(w) and 
hence the orbits of vertices under the function / induce a path-like structure. 

Definition 3 Let (G, /, O) be a geometry. A collection C of (possibly trivial) directed 
paths in G is a path cover o/(G, /, O) if 

(i) . each vertex v S V is contained in exactly one path; 

(ii) . each path in C is either disjoint from I , or intersects I only at its start point; 
(Hi), each path in C intersects O only at its end point. 

A path cover of (G, /, O) is meant to capture the way that information flows within 
the geometry during the execution of a pattern. When |/| = \0\, a path cover is just a 
maximum-size collection of directed paths from I to O which do not intersect (i.e. a 
maximum family of vertex-disjoint I- O paths) which also happens to cover every vertex 
of the graph. 

Every flow function / generates a path cover from its' orbits, and conversely we can 
obtain a natural function / from any path cover. However, not all path covers correspond 
to flows, because the function / from an arbitrary path cover may not have a compatible 
partial order satisfying (Fii) and (Fiii). For example, consider the geometry given by 
alternating input and output vertices on a cycle of length 6: 



- V = {0, 1, 2, 3,4, 5}, I = {0, 2, 4} , and O = {1, 3, 5} ; 

- there is an edge between i, j 6 V iff j — i = ±1 (mod 6) . 

The only two path covers for (G, /, O) are generated by the functions f(x) = x + 1 
(mod 6) and f(x) = x — 1 (mod 6) . In the former case, we have C2 ~ /(co) > C4 ~ 
/(C2) , and Co ~ /(C4) : then, from condition (F/u), we would require Co =^ C2 =3! C4 =<! 
Co . Similarly, in the latter case, we would have Co >p C2 >p C4 fc= Co . However, neither 
of these are possible if =3! is a partial order. 

In order to characterize when a path cover corresponds to a flow, we must introduce 
a condition which prohibits such a cycle of relations between distinct vertices. We do 
this by capturing the conditions which make a path cover give rise to these relations in 
the first place. 

Definition 4 Let T be a family of directed paths in a graph G . A closed walk C = 
cqC\ ■ ■ ■ c m _iCo (indexed over Z m J in G is a vicious circuit for T if for all j € Z m , 

c j-i 7^ Cj+l > and at least one ofcj-i — > Cj or Cj — * c^_i is an arc of T . 

That is: a vicious circuit is a walk that only visits vertices covered by T , never doubles 
back on an edge which has just been traversed, and which traverses an arc of T on 
at least one out of every two consecutive edges. The idea behind considering such a 
structure is that such a vicious circuit can be cut into segments, each consisting of one 
or two edges: 

- One edge segment: an edge xy , where x — > y is an arc of T ; 

- Two edge segment: a path xyz , where x — > y is an arc of T , but y — > z is not. 

The importance of these segments is as follows. Let C be a path cover, and / be the 
injective function naturally coming from it: then x — > y is an arc of C iff y = f(x) . For 
any x and y which are the initial and final points of such a segment, any binary relation 
^ satisfying (Fii) and (Fm) will have x =4 y ■ As vicious circuits can be decomposed 
into such segments, a vicious circuit will then induce the cyclic relations we wish to 
prohibit. This motivates the following definition: 

Definition 5 A causal path cover of (G, I, O) is a path cover which does not have any 
vicious circuits in G . 

The definition above is intended to simply exclude those path covers which induce 
cyclic relations of the sort above. In fact, this definition characterizes flows: 

Theorem 3 (|38|, Corollary 11) A geometry has a flow iff it has a causal path cover. 

6.2 Uniqueness of Causal Path Covers among Families of Disjoint i- O paths 

Theorem|3]tells us that in order to find a flow function / for a geometry (G, I, O) , we 
may replace the condition that / has a compatible partial order satisfying (Fii) and (FhV) 
with the condition that the path cover induced by / lacks vicious circuits. In the special 
case where |/| = \0\ (e.g. for a geometry from the pattern for a unitary transformation), 
we have the following result: 



Theorem 4 (|38|, Theorem 16) Let (G, I, O) be a geometry such that \I\ — \0\ . If 
(G,I,0) has a causal path cover C, then it is the only maximum family of vertex- 
disjoint 1-0 paths in G . 

The reason for this is that any other maximum family of vertex-disjoint 1—0 paths 
either fails to cover all vertices of G, or traverses a different set of edges in G . From 
this, one can prove that a vicious circuit must exist for C . Thus, if C is causal, no other 
maximum family of vertex-disjoint 1-0 paths may exist. As a corollary, because a 
causal path cover for (G, /, O) is unique when it exists, a geometry (G, /, O) has at 
most one possible flow function. 

Using the above result, if (G, /, O) has a flow, we may find it by first finding any 
maximum family of vertex-disjoint I-O paths, which has a known reduction to the 
Maximum Integral Flow problem. 

6.3 Reduction to Maximum Integral Flow 

Informally, an integral flow network consists of a directed graph N, together with a 
single source vertex r and a single sink vertex s , and where each arc x — > y in TV has 
an associated capacity c(x — > y) G N . We are interested in sending come measurable 
commodity from the source to the sink, under the constraints that 

(z). the commodity can only be sent between vertices in the direction of an arc; 
(ii). any amount of the commodity which enters a vertex other than the sink must also 
leave that vertex; 

(Hi), the amount of commodity sent over an arc is at most the capacity of the arc. 

In particular, we are interested in possible network-flow functions u, where u(a) £ N 
gives the amount of the commodity which we send across the arc a . Define 

F u {x) = 

then, the conditions above state that F u (x) = for x ^ {r, s} , and that for each arc 
a , we have u(a) ^ c(a) . The Maximum Integral Flow problem is to find a network- 
flow u such that F u (s) is maximal: this has several known efficient solutions which 
are commonly implemented in programming libraries. We will not discuss them here 
except for the purpose of run-time analysis: interested readers may refer to 1 39 1 or 1 40 1 . 

To reduce finding a maximum family of vertex-disjoint 1-0 paths to Maximum 
Integral Flow, we use the following standard network construction: 

Definition 6 Let a geometry (G, /, O) be fixed. Define a digraph N with the vertices 

- A = {A v | v e I U O c } (out-flow vertices); 

- B = {B v \ v e O U I c } (in-flow vertices), such that A v =B v iffvelnO; 

- a source vertex r and a sink vertex s, 



E 

{w — >x) 



u(w 



(4) 



and with the following arcs: 



- For i £ I, Ai has only the incoming arc r — > Ai ; all other in-take vertices A v have 
only the incoming arc B v — > A v . 

- For uj £ O, has only the outgoing arc — > s ; all other out-take vertices B v 
have only the outgoing arc B v — » A v . 

- If v £ O c and w £ I c , there is an arc A v — > B w iff vw £ G . 

We call N the max-flow digraph of (G, I, O) . 

The purpose of defining this network is to force any network flow function u on the arcs 
of N to describe a collection of paths from r to s in D , from which we can obtain a 
family of vertex-disjoint 1-0 paths in G . We will now describe how this construction 
can be used to find such a maximum family of paths. 

For the case where / and O are disjoint, consider the network that would be obtained 
by adding source and sink vertices r and s to V , replacing each edge xy £ E with arcs 
x — > y and y — > x , adding arcs r — > i and lj — > s for all i £ I and u> £ O . If 
we set the capacities of all arcs in N equal to 1 , then a maximum integral flow in N 
could be represented by a set S of arcs a in N for which u(a) = 1 . From the condition 
F u (x) = for x ^ {r, s} , each vertex other than r or s would have the same number of 
arcs in S which enter it as it does that leaves it. Then, the arcs in S describe a collection 
of edge-disjoint paths from r to s , and this collection of paths is at a maximum when u 
is a maximum integral flow. 

To force verfex-disjointness, we want at most one arc of S to enter each vertex, and 
at most one arc of S to leave each vertex. We can do this by replacing each vertex x by 
two vertices B x — > A x , where all arcs which left x now leave A x , and all arcs which 
entered x now enter B x . Because only one arc leaves B x , at most one arc of S can enter 
B x ; and similarly with arcs leaving A x . Having a maximum family of vertex-disjoint 
paths in N from r to s , we can recover a maximum family of vertex-disjoint paths in 
G by removing r and s , and contracting all pairs {B X ,A X } to a single vertex x £ V . 

For the case where I and O are not disjoint, we can let A x and B x be equal when 
x £ ID O , as long as we do not connect x to any vertices other than r or s . The 
network that results will be the max-flow digraph N described in Definition[6]above. 

Given a geometry (G, /, O) the max-flow digraph for (G, I, O) can easily be con- 
structed in 0(m) time. Then, by the above construction, we can attempt to find a path 
cover for (G, I, O) as follows: 

Algorithm 3 Input: the max-flow digraph N of a geometry (G, I, O) , with arc capac- 
ities set to lfor each arc of N . 

Output: either (i) a path cover C for (G, /, O) and a function f : O c — > I c whose 
orbits induce the paths in C , or (ii) no flow exists for (G, J, O) . 

1. Set u to be a maximum integral network flow for N . 

2. Find a family C of directed paths in G which start at input vertices, end at output 
vertices, and only traverses edges xy such that u is non-zero on the arc A x — > B y . 

3. For each arc A x — > B y of C , set f(x) — y . 

— (i) If the paths ofC cover all vertices of G , return (C, /) . 

— (ii) Otherwise, no flow exists for (G, /, O). 



Whether or not (G, I, O) has a flow, the collection of paths C which is constructed 
will be a maximum family of vertex-disjoint I—O paths; we can check if they cover 
V (and construct the function /) by simply traversing the paths. By Theorem |4j if C 
doesn't cover all vertices, then (G, /, O) has no causal path cover, in which case it also 
has no flow. 

Let n = \V\ , m = \E\ , and k = \I\ = \0\ . The Maximum Integral Flow problem 
on N can be solved in time O(km) by the Ford-Fulkerson algorithm (see, e.g. |39|). 
The time required to construct C and / , and to verify that C is a path cover, will be 
0(n) . Thus the running time of Algorithm|3]is 0(km) when G is a connected graph. 

6.4 Obtaining the Dependency Order of / 

The function / returned as part of Algorithm [3] is necessary, but not sufficient, for 
(G, I, O) to have a flow. We must also discover if / has a compatible partial order 
which satisfies (Fa) and (Ftii). Suppose that / has such a compatible partial order, and 
consider in particular the dependency order =<! of / . Because it is the coarsest compat- 
ible order, =<! can be defined as the transitive closure of the binary relation d on V , 
given by 

xdy [y = x] V [y = f(x)\ V [y ~ f(x)} , (5) 

which we will call the dependency relation: we will say that y depends on x if x d y . 
We can then reduce finding the dependency order to the Transitive Closure problem, 
which is concerned precisely with finding the transitive closure of binary relations. 2 If 
there is a cycle of dependency relations, then the transitive closure of d will not be a 
partial order, in which case / is not a flow function; otherwise, / is a flow function, and 
the transitive closure will be the dependency relation of / . 

Algorithmic is a Tarjan-like algorithm for detecting cycles in the dependency re- 
lation d (given the function / and graph G which determines d ), and computing the 
transitive closure if no such cycles exist. It also takes the path cover C as input. 

Algorithm 4 Input: A geometry (G, I, O) , and the output (C, f) of Algorithm^ 
Output: either (i) the dependency order =^ of f, or (ii) no flow exists for (G, I, O) . 

— Let sup be an array mapping pairs (x,P) G V x C to the supremum (i.e. the least 
upper bound) of the set {x} in the path P , or the empty set if there is no least upper 
bound. 1 Initialize svp[x, P] = for all (x, P) G V x C . 

— Let status be an array mapping x £ V to one of {none , pending, fixed}; 
where statusja;] = none indicates that no supremum of {x} has yet been dis- 
covered, and status [x] = fixed indicates that all suprema of {x} have been 
discovered. Initialize status [x] = none for all x G V . 

2 The Transitive Closure problem is usually expressed in terms of directed graphs, but doing so 
doesn't provide any additional insight into the problem, and is at any rate equivalent to the 
problem for binary relations. 

3 In the literature for the Transitive Closure problem, this is known as a chain decomposition of 
the directed graph represented by the relation d . 



- For each vertex x E V with status [x] = none perform the following recursive 
algorithm to set the suprema of {x} : 

1. Set status [x] = pending. 

2. For the path P x which contains x , set sup[x, P x ] = x . 

3. For each vertex y which depends on x (aside from x itself): 

• status [y] — none, recursively determine the suprema of {y} . 

• if after this status[y] 7^ fixed, (G, /, O) has no flow; abort. 

• otherwise, set sup[x. P] — min( sup[.r. P] , sup[j/, P] ) for each path 
P € C (taking the minimum with respect to distance from the start of P). 

4. Set status [x] = fixed. 

- Once status [a;] = fixed/or all x E V , return the array sup. 

This algorithm determines if there is a cycle of dependencies by performing a depth- 
first search, and checking if it traverses a sequence of relations xdyd---dx'm the 
course of finding the elements which are greater than x in the dependency order. If it 
does, it will discover that status [x] = pending on the second time it visits x : then, 
it may as well abort, as / cannot be a flow function and no dependency order exists. 

The array sup returned as output characterizes the dependency order =^ : we have 
x =4 V if and only if sup [a;, P y ] ^ y , where P y is the path containing the vertex y , 
and the inequality ^ is with respect to the distance from the initial vertex of P y . We 
can compute this easily by constructing functions which maps each vertex y G V to 
an index for P y , and another function mapping each vertex x E V to the distance of 
x from the initial point of P x . This can be considered as a part of building the path 
cover C in Algorithmic (indeed, such functions characterize the path cover C and may 
be considered to be how C is represented). 

Let n = \V\ , m = \E\ , and k = \I\ = \0\ = \C\ . The time required to fix the 
values of sup[x, P] for each path P E C is 0(k deg(x)) . Summing over all x E V , the 
cost of recursively obtaining the suprema for all of V is O(km) . This is then the cost 
of Algorithm [4] as a whole, as initializing sup and status takes a total of O(kn) C 
O(km) time for G connected. 

7 Discussion and Conclusion 

We have presented a decomposition of unitary maps into three successive operations 
P, <!>, and R. The first one represents the familiar preparation map and expands the 
computation space by introducing auxiliary qubits; the second one is diagonal in the 
computational basis and has only unit coefficients; and the last is a restriction map that 
contracts back the computational space into the chosen output space. It is important to 
emphasize that both R and P have a very simple structure, and hence the decomposition 
suggests that the whole quantum computing part of an algorithm is encoded in the phase 
map operator. 

The restriction map R does not correspond directly to a meaningful physical trans- 
formation, so this decomposition is a mathematical artifact. However, we have illus- 
trated how one may attempt to find a one-way patterns which effectively implements 
such phase decompositions. We have also shown how one may attempt to find a phase 



map decomposition directly. This constitutes a method for the implementation of uni- 
tary maps in the one-way model which completely bypasses the circuit model. It is 
hoped that the phase map decomposition may produce more efficient implementations 
in particular instances. 

An interesting restriction which may be imposed is to restrict the choice of the 
allowed angles. For example, we might restrict to multiples of 7r/4, in the hope of 
staying within the world of Pauli measurements and the so-called 7r/4 magic prepara- 
tions, which may help with implementations 1171 . In this restricted case, the Equations 
of Lemma|2]may no longer have an exact solution in general; however, one does not 
expect an exact solution in this case, since the obtained model is only approximately 
universal. It would be interesting to see whether one can compute down an approximate 
phase map decomposition such that U « R o <P o P with arbitrary precision. 

Another question is whether there is a reasonable relaxation of the definition of a 
pattern, by allowing new commands which would allow additional possible phase maps. 
How much a non-trivial extension would complicate the theory (e.g., standardization, 
flow) and the implementation (e.g., requiring only one simple primitive for entangle- 
ment) of one-way patterns is not clear. 
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